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Abstract 

For the phenomenon of Bose-Einstein condensation we construct the 
canonical pair of field operators of the Goldstone Bosons explicitly as 
fluctuation operators in the ground state. We consider the imperfect 
Bose gas as well as the weakly interacting Bose gas. We prove that a 
canonical pair of fluctuation operators is always related to the order 
parameter and the generator of the broken symmetry fluctuations. 
We find that although the first one has an anomalous behaviour, the 
second one is squeezed by the same inverse rate. Furthermore, we 
prove that this canonical pair separates from the other variables of 
the system and that it behaves dynamically as oscillator variables. 
Finally the long wavelength behaviour of the spectrum determines 
the lifetime of this pair. 
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1 Introduction 



The phenomenon of spontaneous symmetry breaking (SSB) is a representa- 
tive tool for the explanation of many phenomena in modern physics of field 
theory and statistical mechanics of many-body theory. The analysis of SSB 
goes back to the Goldstone Theorem which has been the subject of much 
analysis. It is proved that for short range interactions in many-body systems 
SSB implies the absence of an energy gap in the excitation spectrum || §]. 

For long range interactions the SSB has also been studied extensively. In the 
physics literature the phenomenon is known as the occurence of oscillations 
with frequency spectrum taking a finite value u ^ at k = j|, |], |(J . Differ- 
ent approximation methods, typical here is the random phase approximation, 
yield the exact computation of these frequencies. For some mean field mod- 
els, the BCS-model f7|, the Overhauser model f|, the anharmonic crystal 



model ||, and for the Jellium model [IU|, we were able to give the math- 
ematical status of these frequencies as elements of the spectrum of typical 
fluctuation operators (see jD], |T2]|). 

The typical operators entering in the discusssion are the generator of the 
broken symmetry and the order parameter. In physical terms expressed, it 
is the charge or density operator and the current operator. Their fluctu- 
ation operators form a quantum canonical pair, which decouples from the 
other degrees of freedom of the system. As fluctuation operators are col- 
lective operators, they describe the collective mode accompanying the SSB 
phenomenon. Hence for long range interacting systems, we realised mathe- 
matically rigorously in these models, the so-called Anderson theorem |fTB| , [H|j 
of 'restauration of symmetry', stating that there exists a spectrum of col- 
lective modes uo{k — > 0) 7^ and that the mode in the limit k — > is the 
operator which connects the set of degenerate temperature states, i.e. 'ro- 
tates' one ergodic state into an other. We conjecture that our results of 
J7|, |8], 10] can be proved for general long range two-body interacting systems 
as a universal theorem. Anderson did formulate his theorem in the context 
of the Goldstone theorem for short range interacting systems, i.e. in the case 
u{k — ► 0) = of absence of an energy gap in the ground state. 

Of course one knows that there is no one-to-one relation between short range 
interactions and the absence of an energy gap for symmetry breaking systems 
(see e.g. ||). The imperfect Bose gas is an example of a long range interacting 
system showing SSB, but without energy gap. In this paper we realise the 
above described program of construction of the collective modes operators 
of condensate density and condensate current, as normal modes dynamically 
independent from the other degrees of freedom of the system. We consider 
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the whole temperature range, the ground state included. 

In particular the ground state situation is interesting, because it yields a non- 
trivial quantum mechanical canonical pair of conjugate operators, giving an 
explicit representation of the field variables of the socalled Goldstone boson. 
One can consider this result as a formal step forward beyond the known 
analysis of the Goldstone phenomenon. 

Moreover in Section |], we extend this result tot the weakly interacting Bose 
gas of superfluidity. It is interesting to remark that here the situation is 
intrinsically different in the sense that only the condenste density mode is 
spontaneously broken. One checks explicitly that the density fluctuation 
operator and the order parameter fluctuation operator do not form a non- 
trivial pair, but the condensate density and the order parameter fluctuation 
operators do. 

Hence in both models the fluctuation operators of the generator of the broken 
symmetry and of the order parameter form a non-trivial canonical pair. The 
latter one shows off-diagonal long range order, therefore the density-density 
correlation can not share this property. This can be interpreted as that 
a spontaneously broken symmetry behaves like an approximate symmetry. 
The explicit construction of the canonical pair amounts to the realisation of 



'restauration of symmetry', an idea put forward by Anderson |T^, [14]. 

Furthermore, for both models, we prove that the canonical pair of Gold- 
stone fluctuation modes separates dynamically from the other variables of 
the system and behaves like harmonic oscillator modes with a frequency pro- 
portional to the condensate density, i.e. this phenomenon disappears if no 
condensation is present. It turns out that this pair of variables has a lifetime 
in the long wavelength limit which is determined by the long wavelength 
behaviour of the spectrum of the system. 



2 Fluctuation operators 

We want to study different models of a Bose gas in which there is breaking 
of the gauge symmetry. In general, a system of identical bosons of mass m 
in a cubic box A C M" of volume V = If , v > 3 with periodic boundary con- 
ditions for the wave functions, is described by the full two-body interaction 
Hamiltonian 

H L = Y, e kO* Ltk a LM + — v (<l) a L,k+<i a L,k>- q a L,k>aL,k ~ HlN l , (1) 

k ZV q,k,k' 
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where the sum runs over the set A* = and = a L k are the boson 
creation/annihilation operators in the one-particle state vpL,k(x) = V~ l ^ 2 e tk ' x , 
x e A, k G A*, i.e. 

aL,fc = y a(x)^ 7T y 5 -cia;, (2) 

[a(x),a*(y)] =8{x-y) 

and = f R „ e~ iq ' x (f>(x)dx, (ft is the periodically extended two-body inter- 
action potential. 

The generator of the gauge symmetry is the total number operator Nl, with 
generator density a*(x)a(x): 

Nl = / a*(x)a(x)dx. 

J A 

The common choice of order parameter is V~ 1 ^ 2 a L0 , or taking a self-adjoint 
combination 

i if 

so the order parameter density is given by A^{a*{x) — a(x)). One has of 
course. 

[N L ,A L ] = -=(a* Lfl + a Lfi ) v -^ zy^cosa, 

V — V 

where p is the density of the condensate and a the phase, i.e. 

We are here interested in the behaviour of the g-mode fluctuation (q ^ 0) 



|T5[ of this generator and order parameter, i.e. 



F L , q (N) = ^ I a U)ai.r)< ■"-<!,■ (3) 



A 



F L , q (A) = -^-j-JJa*( x )- a (x))e«*dx, (4) 
which satisfy the same commutation relation as Nl and Al'- 

[F L , q (N),F L ,- q (A)} = -±=(al fi + a L , ) V -=^ zy^cosa. (5) 
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In fact we take a sequence 7^ qi £ A* converging to q so that there is no 
need to substract expectation values, since for g G A*, J A e %q - x dx = V5 qt o, and 
so that we can also write 

F L,q(N) = T^Y, a L,k+ q a L,k (6) 

k 

F L , q {A) = * (o^-a^-,). (7) 



Our first goal will be to define these operators in the thermodynamic limit 
L — > 00. This will be done via a central limit theorem, as defined in [0], |T2| |. 
Afterwards we will be interested in the long wavelength - low frequency limit 
q — > in which collective behaviour is to be expected. In this limit q — > 
we give a connection with the abstractly studied fluctuation operators [16] 
of the type 

F s (0) = lim F L , 5 (0) = -4-r I (0(x)- < 0(x) >)dx, (8) 

where O is some operator density and 5 a critical exponent describing the 
degree of abnormality of the fluctuations of O, defined by existence of the 
variance. If 5 > there is ODLRO, if 5 < 0, the fluctuation is squeezed. 

In an interacting Bose gas this q — > behaviour will be mainly determined 
by the spectrum E q of the Hamiltonian. The density fluctuation F L ^ q (N) 
has another very important property, its commutator with the two-body 
interaction part of the Hamiltonian vanishes: 

[tf L ,F L ,,(JV)] = 0, (9) 

where 

U L = v ((l) a L,k+ q a L,k'~ q a L,k'a L , k . 

ZV q,k,k> 

This is easily seen as follows. Commute in Ul, d* L k ,_ q with ai,k- This gives 

1 1 
U l = 7f{7 zl v ( ( l) a L,k+ q a L,kal tk/ _ q a Lr - 7 p(7Y, v ( ( l) a L,k+ q a L,k+ q - 

ZV q,k,k' ZV q,k 

In the first term, separate the term q = from the rest, use translation 
invariance in the second term, and observe that Ul can then be written as 

Ul = \ Y,v(q)FL, q (N)F L ,- q (N) + ^iVf - ^(0)iV L . (10) 
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From this expression, ([)]) is obvious. 

In physics, one encounters essentially two types of Bose condensed systems, 
namely those with a quadratic excitation spectrum (E q oc \q\ 2 , \q\ small) and 
those with a superfluid, linear spectrum (E q oc |g|, \q\ small). We will treat 
in detail an example of each of these cases. 



3 The Imperfect Bose Gas 

3.1 The model and equilibrium states 

To make things more concrete, we consider as a first example the imperfect 



or mean field Bose gas [17, |TB|, specified by the local Hamiltonian Hi with 
periodic boundary conditions JjlJ: 

H l = T l -hlN l + ±NI (11) 

where A G M + , A the centered cubic box of side length L in R", v > 3, 
A* = ^j- r L v '. Remark that, apart from a shift in the chemical potential, this 
Hamiltonian can be obtained from ([]]) by taking v(q) = if q ^ in fllPp. 

Talking about the thermodynamic limit, we mean L — > oo under the con- 
straint that for all L 

(12) 

where p is any positive number standing for the average density of particles, 
ujl is the canonical Gibbs state for ([TT|) at some inverse temperature j3. It is 
proved that u)p(.) = limxo;i(.) exists as a space homogeneous state on 
the algebra of polynomials in the creation and annihilation operators. It is 
proved that there exists condensation in the zero (k = 0) mode state if p is 
large enough and T is small enough. 

The phase transition is accompanied by a spontaneous breaking of the gauge 
symmetry of flTT| ) in the sense that 

hm u L (^a(f) \ = p /(0) , po > 0, (13) 



where a(f) = J f(x)a(x)dx for / e C^°(W), f is the Fourier transform of 
/ and po is the condensate density. It is proved in |20| that (|13"D implies, 



amongst other things, the breaking of the gauge symmetry. 
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The limit chemical potential is given by 

// = Xp, 



(14) 



and the dynamics coincides with the dynamics of the free Bose gas. The 
limit Gibbs state has the following form: for all local observables A, 



with 



and 



u a ( e i(a(/)+«-(/))) =exp 



--(/,#/) + 2v^|/(0)|cosa 



(iff) (*) = icoth(^)/». 



(15) 

(16) 
(17) 



The states u/g(a G [0, 2ir]) are the extremal equilibrium state components of 
Up with the property that 



( a L \ 



Poe 



(18) 



and as operators in the GNS-representation of lu^, one has also 



L L,0 



lim 

z^oo yi/2 



Poe 



(19) 



Remark also that the states ujp(a G [0, 2ir]) are quasi- free states, making the 
computation of expectation values straightforward. 



3.2 Collective Goldstone modes 

We now turn our attention to the density and order parameter fluctuations. 
We consider our system to be in one of the extremal equilibrium states u/g, 
for some a G [0, 2ir], and without loss of generality, we take a — 0, and 
denote this state again by up. 

For notational convenience, if po 7^ 0, denote 

PL ' q = = Fl ^- 



Then we have 



[pL, q , A £ = 2v /^y K,o + «l,o), 



7 



and by ([IJ): 
More generally 



lim \p L A L ,_ q ] = i. (20) 

L—>oo 



lim \p L q , A L g,] = iS q q/ . (21) 

L — >oo 



Let us first calculate the variances of pL, q and Al^. 
Proposition 1 We have for q, q' ^ 
(i) 

lim U)p(p L ,qPL,-q>) 
L— >oo 



1 j _ 1 /3e g 1 f dk 



(ii) 



Sao' \ - COtll 1 , ; — 

qA \ 2 2 2p A" (2vr)^ e^ e *+« - 1 1 - 



lim ^(Ai^Ai-g/) = ^g/icoth^. 

L— >oo Z Z 



Proof. The proof is a straightforward calculation using the quasi-freeness of 
the state ujp, e.g. 

Uf3(pL, q pL,-q) 

= o — T7zC^( a I,fc+g a ^,fc+9) u; /3( a i,fc a L.fc) 
^PoV k 

= 7^7 (uf3(a* L:q a Liq )uj(3(a Lfi a* Lfi ) + ujp(a* Li0 a LjQ )ujp(a L - q a 



L,—q, 



In the limit, this becomes 



lim U/3(p L , q PL-q) = o coth + 



1 xl (3e q , 1 [ dk 1 1 



l-oo p^>«^> 97 2 2 2p y R « (2 7 r) l/ e^+« - 1 1 - e~^' 

The other case is even easier. 



□ 



From this already a few conclusions can be drawn. First of all, consider the 
integral in the most relevant case v = 3: 

1 r dk 1 1 
2p^ Jr3 (2vr) 3 e^+i - 1 1 - e-P'* ' 

Letting q — > 0, this integral clearly diverges due to the contribution of the 
neighbourhood of k = 0. Near k = we can write it (up to constants) like 

r dk 1 
J (2vr) 3 (k + qfk 2 ' 

Taking q e.g. along the z-axis and changing the variable k to k' = A, it 

can be seen that this integral diverges like Iq]' 1 - Since coth ^ diverges as 
\q\~ 2 for q —>■ 0, we see that for small q the variance of pi, q is completely 
dominated by the coth-term. 

This divergence implies that we should renormalise both pi iQ and Ai. q in 
order to get a nontrivial limit q — > of their variances, i.e. 

Ph,q —> PL,q = \q\PL,q, 

But this implies that the commutator 

lim [p L ,q,^L,-q] = i\q\ 2 , 
vanishes in the limit q — > 0. 

On the other hand, if one considers the ground state situation (limit (3 — > oo), 
the q — > analysis yields that the variances 

timUoo(pL,qPL,- q ) 

and 

are both finite, and the commutation relation between p^ A and A^ _ q is non- 
trivial and canonical. 

This is not surprising, as one expects true quantum effects on the level of 
fluctuations only in the ground state. Critical quantum effects are hidden 
behind the temperature (T > 0) fluctuations. 

Moreover, it is tempting to identify this renormalisation in q with the expo- 
nent 5 of (§) via the relation |g| oc L^ 1 . This relation of course being given 
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by the fact that the first non-zero g-level in finite volume is \q\ = 2nL~ l . 
In that case we obtain for the density fluctuations F^^N) that 5 = 1/3 in 
the condensed phase. In the normal phase, the coth-term would be absent 
and the integral would be convergent also for q = because e^ tk would be 
replaced by e^ ek ~ a \ with a < 0, hence 5 = 0. At the critical point, the 
coth-term would still be absent but the integral would now be divergent like 
as shown before. This would then give 5 = 1/6. These three values for 
5 are exactly the ones calculated in 



Since we will be interested in quantum effects on the level of macroscopic 
fluctuations we will restrict ourself from now on to the ground state (from 
now on denoted uj). We redefine pL, q and Al, q as self-adjoint operators, and 
we make the (arbitrary) choice of taking the cos- fluctuation: 



pLq = —== I a*(x)a(x) cos(q.x)dx (22) 
y/poV J A 



Ai, q = —j= \ (a*(x) — a(x)) cos(q.x)dx, (23) 



or in momentum space 



PL, q = 2 Jp~oV ^ a * L *+i aL ' k + a *L,k~q a L,k) ( 24 ) 



k 



A L, q = ^[ a *L, q + o* L _ q - (a L>q + a L - q )\, (25) 
where the normalization is chosen such that 

lim \p L q , A L ql ) = i5 qtq >. (26) 

It is easy to check that with these definitions 

hm^pij = \ (27) 

= \- (28) 

The rest of this section is devoted to the more mathematical aspects of the 
realisation of the different fluctuation operators as central limits of operators. 
The less mathematics minded reader can skip this part at a first reading and 
proceed immediately to section |3~5. 



Let T be the family of complex continuous functions f(k, k') of two variables 
k, k' E W, satisfying 

f(±k,±k') = f(k,k') (29) 
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and 

7(M0 = /(*',*). ( 3 °) 

With later applications in mind, define for f,g e J 7 , pL,q(f) and AL, q (g) by 



PL,q(f) 



+f(k-q,k)at k _ q a Lik 



(31) 



S(9> °) ( a L,g + a I,-J - ?) K,9 + a L,-q) ■ (32) 



Condition (^) ensures the self-adjointness of these operators. Then define 
operators F L , q (f,g) by 



F L , q (f>9)=PL,q(f)+ A LM' 



(33) 



Proposition 2 For f,g G 

hm W (F L , g (/^) 2 ) = i|/(g,0)+^(g,0)| 2 . 



(34) 



Proof. This is a simple calculation using the quasi-freeness of the state uj. 



□ 



In the GNS-representation (H w , tt^, Q w ) of the state uj a scalar product is 
defined by 

and the associated norm is denoted by ||.||a,. 

Denoting 7r u (F Lt g(f, g)) again by F L , q (f,g), we then have: 

Proposition 3 (A BCH-formula) Let fi,gi 6 J 7 , i = 1,2, i/ien 



lim 



e iF L , q (h, gi ) e iF L Jf 2 ,g 2 ) _ e i{F L , q {h, gi )+F L , q {h m )) e -\[F L , q {f im ),F L , q U 2 , g2 )] 



(35) 



Proof. See Appendix [A]. 



□ 
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This result should be compared with the Baker-Campbell-Haussdorff formula 
which states that for two operators A, B whose commutator is a complex 
number: 

e A e B = e ^ A+B) e~^ [A ' B] . 

This proposition tells us that in a weak sense this BCH-formula remains true 
for our fluctuation operators, whose commutator becomes a complex number 
in the thermodynamic limit. Since we are studying fluctuations of unbounded 
operators, the BCH-formula is only true in the GNS-representation of uj. For 
fluctuations of bounded operators, the BCH-formula holds in a much stronger 
sense, independent of the state (see |T5|| ). 

On the complex vectorspace V of complex linear combinations of elements 
from (JF, J 7 ), define a sesquilinear form < -|- > q by 

< fi, 9i\f2,92 > q = lim Lo(F Lq (f 1 ,g 1 )*F Lq (f 2 ,g 2 )) 

= \(fi(q, 0) + i gi (q, 0)) (f 2 (q, 0) + ig 2 (q, 0)) , (36) 

and extension to the whole of V by linearity. 

This form is positive and satisfies the Cauchy-Schwarz inequality by the 
positivity of uj and the Cauchy-Schwarz inequality for the state to. 

Separating the real and the imaginary part of the restriction of < -|- > q to 
the real subspace (J 7 , JF) of V, i.e. 

i 

< fi,9i\f2,92 > q = s q (f 1 ,gi\f 2 ,92) + -<r q (fi,gi\f2,g 2 ) , (37) 

defines a real bilinear positive symmetric form s q and a symplectic form o q . 
(A form a is called symplectic if cr(x,y) = —a(y,x).) 

The symplectic form a q satisfies 

lim [F Liq (f 1 ,g 1 ),F Ljq (f2,g 2 )] = i° q , (38) 

L—>oc 

where the limit is taken in the GNS-representation of to. 

The following proposition is the crucial Central Limit Theorem for the oper- 
ators F L>q (f,g). 

Proposition 4 (Central Limit Theorem) For f,gGJ r ,tEM., 

lim u ( e itFL ^ f ' 9) ) = e -#^(/,9l/,s). (39) 
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Proof. Although a similar theorem could also be proven for temperature 
states, we will only do it for the ground state, since that is really all we 
need. In that case, using the quasi-freeness of the state and the fact that all 
particles are condensed into the zero-energy state simplifies the proof. The 
details can be found in Appendix [B]. 



□ 



The C*-algebra of the canonical commutation relations over (H,o~), with H 
a real linear space and a a symplectic form, written as CCR(H,a), is by 
definition a C*-algebra generated by elements {W(f) : / G H} such that 



(i) w(-f) = w(fy 

(ii) W{f)W(g) = el"V,9)w(f + g). 

Condition (ii) tells us that W(f)W(0) = W(0)W(f) = W(f). Hence W(0) 
is the unit of the algebra and it follows that W(f) is a unitary for every /. 
For an elaborate discussion of the CCR, we refer to |p2|. 



Proposition 5 (Reconstruction Theorem) The linear functional 

Cd q (W q (f,g)) = e -K(/.9l/.s) (40) 

defined on the algebra CCR((J r ,J-'),o- q ), is a quasi free state. 
More explicitly, we have for all (fi, gi) G (J 7 , J-), i = 1, . . . , n, 

lim u) ( e «^<A.«) • • • e iF ^ f ^"A = W {W q {h, g x ) ■ ■ ■ W q (f n , g n )) . (41) 

The state uj q is regular and hence for every (f, g) there exists a self-adjoint 
Bo sonic field <& q (f,g) in the GNS representation (H^ , vr^q , Q^i ) such that 

^(W q (f,g)) = e^^. (42) 

This implies that in the sense of the central limit the local fluctuations 
converge to the Bosonic fields associated with CCR((J r ,J r ),a q ): 

CUT - hm F L , q (f,9) = * 9 (f>9)- (43) 

L — ►oo 
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Proof. See Appendix |CJ 



□ 



The following definitions now clearly make sense: 

p q = $,(1,0) = CLT- limp X)? , (44) 
A q = $ g (0, 1) = CLT - Km A L>q . (45) 

In the same spirit as this central limit, we now define a limit q — > of the 
operators $ q (f,g). 

Define a sesquilinear form < -|- > on V by 

< fi,9i\h,92 >= I™ < fi,gx\f2,92 > q , (46) 

q-+0 

and a real linear form s and a symplectic form a in the obvious way: 

gi\hi92) = (fx, 9i\f 2 , 92) (47) 

cr(fi, 9i\f2,9z) = Km<T q (fi,gi\fz,g 2 ). (48) 

We then get the limit (q — > 0) result: 

Proposition 6 (Reconstruction Theorem 2) The linear functional 

u>(W(f,g)) = limu« (W q (f,g)) = e "H/^) (49) 

defined on the algebra CCR^J 7 , J 7 ), a), is a quasi free state. 
More explicitly, we have for all (fi, g,i) G (CF, !F),i — 1, . . . , n, 

lim lim u ( e iF ^ fl > 9 ^ ■ ■ ■ e iF ^ fn ' 9n) ) = Co (W(fi, g x ) ■ - ■ W(f n , g n )) . (50) 

q— >0 L— >oo V / 

The state Co is regular and hence for every (f, g) there exists a self-adjoint 
Bosonic field $(/, g) in the GNS representation (Hz, 7T&, Q&) such that 

^(W(f,g)) = e^\ (51) 



This implies that in this sense of the limit q — > (\5Q), the fluctuations & g 
converge to the Bosonic fields associated with CCR((CF,CF),a): 

$(f,g) = lim $,(/,<?) = CLT - lim lim F L>q (f,g). (52) 

q^O q—fU L— >oo 
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Proof. This is just a matter of taking the limit q — > in the different steps 
of the proof of the previous Proposition. 

□ 

Specifying again to our original operators: 

p = lim p q = $(1, 0) = CLT - lim lim p L Q , (53) 

q^O q^O L^oo 

A = lim Aa = $(0, 1) = CLT - lim lim A L . (54) 

The algebra of macroscopic fluctuations CGRi^T \!F),a) is a coarse grained 
one, i.e. different microscopic observables can have the same macroscopic 
fluctuation operators. To describe this mathematically, introduce an equiv- 
alence relation ~ on V by 

{h,9x) ~ (/2,0O < h~ h,9i- 92\h- h,gx~ 92 >=0. (55) 

Another way of stating the equivalence relation is of course 

Cfi, 9i) ~ (/2, #2) ^* lim lim w (F L>q (fx - f 2 , gi - # 2 ) 2 ) = 0. (56) 

q^O L^oo \ ' 

We then have the following result: 

Proposition 7 For fi,gi & F,i = 1,2, the following are equivalent: 

(1) (fugi) ~ (h, 92) 
(a) $(fi,gi) = $(/2,0O. 

Proof. See Appendix 0. 

□ 

A simple example: take / G J 7 , and define Jf by (Jf)(q,0) = —if(q,0), 
(J/)(0, g) = «/(0, g) and (Jf)(k, k') = for all other values of k en k'. Then 
(/> 0) ~ (0) Jf) or i n other words 

CLT - lim lim p L Jf) = CLT - lim lim A L JJf). (57) 

q— >0 L^oo ' L^oo ' 
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3.3 Dynamics of the collective Goldstone modes 



In this section we will derive a dynamics on the level of the macroscopic 
fluctuations. This dynamics will of course be induced by the microdynamics. 
Therefore we start with calculating 



i[H L ,pL, q 



with e(k, k') 
And also 

i[H L ,A Ljq 



i{tk ~ efc')- 

1 



T/2 ^2 [( C fc+<? _ e k)a*L,k+q a L,k + (ffe-g ~ ^k) 0*L ,k~ q a L ,k 



(58) 



A 



V 



l L,q 



L-q 



+ a L,-q + CLL,, 



^7 ( a L,q + a l,-q - ( a L,~q + a L , q ) 
A L,q(t) - \ [y NL ~ { a L,q + a L,-q + ^L-q + d Ljq ) 



iX 



a L,q + a L-q ~ \ a L-q + 0>L,q, 



(59) 



The second and the third term on the r.h.s. converge to zero as L —>■ oo, even 
as operators in the GNS repesentation of uj, so they are of no importance. 

Remark that both lim^-^ uj(pi tq (e) 2 ) oc e 2 and lim^oo u(AL, q (e) 2 ) oc e 2 , so 
it is natural to define a macroscopic dynamics by 



i[H,p] = CLT — lim lim % 

q— >0 L— >oo 



—Hi, p L q 



= CLT — lim lim p L „ — e 

9-0L-,oo r ' q \ tq J 

= CLT — lim lim At „ 

= A 

where we have used equation (|57D to go from the second line to the third. 
Analogously, 



i\H,A] = CLT -lim lim i 



1 



~H L ,Ar 



-CLT - lim lim A L , 

q — >0 L—>oo ' 
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-CLT — lim lim or a 



again using (|57|). 

Hence we have found a canonical pair of observables p and A, satisfying 

\p,A]=i, (60) 

which dynamically decouple from the other degrees of freedom of the system, 
with the dynamics given by 

i[H,p] = A (61) 
i[H,A] = -p. (62) 

So H is the harmonic oscillator Hamiltonian with frequency 1: 

H = \ (f + A 2 ) • (63) 
The virial theorem is also satisfied, i.e. 

Co (p 2 ) = Co (A 2 ) . (64) 

Remark that to go from the microdynamics Hl to the macrodynamics H we 
had to rescale the Hamiltonian with el" 1 . This should actually be seen as a 
rescaling of time 

~ t 
t^t = — , 

indicating that for small \q\ the typical lifetime of a (density) fluctuation with 
wave length \q\~ l is of the order \q\~ 2 , becoming infinite in the limit q — > 0. 



4 The weakly interacting Bose gas 
4.1 The model and equilibrium states 

Our second model is a model of superfluidity, i.e. with an excitation spectrum 
Eg linear in |g| for small q. Such a model is provided by f23 |. Its Hamiltonian 
is given by (we take v = 3 throughout this section) 

h l{c) = e k a l,kaL,k + - v ( k )( a l,k a L,-k c2 + c 2 a L ^ k a LM ) 

k / fc^O 

+ |c| 2 £ v(k)al tk a Lik + V M N l _ Aii jv £j (65) 
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supplemented with 

c= \imcu L (V- 1/2 a Lfl ), (66) 

L— >oo 

where ujl is the Gibbs state at some inverse temperature j3 corresponding to 



This Hamiltonian is in fact the original Bogoliubov Hamiltonian for a weakly 
interacting Bose gas, with an extra term ^-N 2 that ensures the supersta- 
bility of the model. 

Again we take the thermodynamic limit under the constraint 

lim ^-u L (N L ) = p. 

L— >oo y 



It is proved in [£3[ that there exist solutions up = lim^ooU^, for (3 and p 



large enough, such that 



limw /J (r%) = c^0 ) 



and we will restrict ourself to these solutions. 

To describe these equilibrium states we need a Bogoliubov transformation of 
the operators a\ k into new creation and annihilation operators b\ k : 

aL,k = b L)k cosh a k + b* L _ k sinh a k (67) 

a-L-k = b L _ k cosha k + b* Lk smha k , (68) 

where ^ 

tanh2a fc = ^MtTTTT - ( 69 ) 



The limit Gibbs state has the same form as in the imperfect Bose gas: for 
all local observables A, 



1 r 2n 

with 



MA) = i ^(A)da, (70) 



u ° ( e W)+"*(/))) =exp 



-(/,#'/) + 2z| C /(0)|cosa 



(71) 



and 

(X7)(A:) = icothf^/(fc) l (72) 
with &"(/) the corresponding Bogoliubov transformation of a"(/) and 

£ fc = ^(e k +\c\Mk)) 2 -(\c\Mk)) 2 = ^ k (e k + 2\c\Mk))- (73) 
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This is the famous Bogoliubov spectrum which for small k behaves like 



1/2 

i \k\. 



m 



The states Ug(a G [0, 2tt]) are the quasi- free extremal equilibrium state com- 
ponents of up with the property that 



and as operators in the GNS-representation of one has also 



}SLv$ = \^"- (75) 



For more details we refer to [E3 . 



As in the imperfect Bose gas, we want to study the density and order parame- 
ter fluctuations, Fi iq (N) and Fl a (A). However here, the Hamiltonian Hl(c) 
(|65f) is a truncation of the full Hamiltonian flU), and due to this truncation, 
H L (c) is no longer gauge invariant, i.e. 

[H L (c),N L }^0. (76) 

The invariance which is left is 

[H L (c),N L ,o) = 0, (77) 

with Nl q = a* L cil,o This means that the spontaneously broken symmetry 
accompanying the phase transition from c = to c ^ is not the gauge 
symmetry generated by Nl, but the symmetry generated by N^q. 

One example of the implications of this is the following. In the physics 
literature, the quantity lim^oo < F L!q (N)F L _ q (N) > is known as the static 
structure function, usually denoted S(q). It has been known for a long time, 
both theoretically and experimentally, that at zero temperature this function 
behaves linearly in q for small q: S(q) oc \q\ (see e.g. This linear 



behaviour is essentially due to the fact that [Ul, F L ^ q (N)] = so that 

<[F L>q (N),[H L ,F L) „ q (N)}\ >oc \q\ 
However in our model 

[H L (c),N L ]?0, 
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and very much related to this, 

[U L (c),F Liq (N)}^0, 

and indeed it is easy to calculate that here lim^o S(q) = const 7^ 0. 

Because of (|76|) and (|77|) we expect that this unphysical behaviour is remedied 
when we replace the total density fluctuations F Lq (N) by condensate density 
fluctuations FL tq (N ). However since Nl$ can not be written as the integral 
over some condensate density, it is impossible to define FL tq (N ) as a usual 
fluctuation operator. What we want to show now is that it is possible to 
find a fluctuation operator FL q (N ) which behaves mathematically like one 
expects for a fluctuation operator of the generator of a spontaneously broken 
symmetry (i.e. we will derive a similar structure as in the imperfect Bose 
gas) and moreover gives the correct physical behaviour (like e.g. 

lim < F Lt g(No)F L .- q (N ) >oc \q\ 

L — >oo 

for small q). 

In momentum space we have 

F ^ N ) = ^E a W^- (78) 

k 

This consists of two parts 

F L,q( N ) = -f7T72 ( a L,q a L,0 + a lfl a L,- q ) + £ a l,k+ q a L,k- (79) 

V V k^0,k+ q ^0 

The first part 

{ a L, q a L,0 + a*L,0 a L,-q) 

is the part of Fl^(N) which contains the ground state operators a^ . It 
is clearly the fluctuation of the zero-mode particle density, fluctuating to a 
fixed mode and back to zero. The other part is the fluctuation of the excited 
modes among each other. Therefore it is natural to define 

F L , q (N ) = p-j^ (a*L, q a L fl + a% a L - q ) . (80) 

The truncation of Fi. q (N) to Fl^(Nq) reminds very much the spirit behind 
the truncation which led to the Bogoliubov approximation of the full Hamil- 
tonian. We can even show how closely those two are related. Take the 
interaction part Ul of the full Hamiltonian ([I]) and write it as in ([T0|): 

Ul = \Y. v{k)F L , k {N)FL t -k{N) + V -^Nt - U{0)Nl. 
z k^o AV A 



20 



Truncate this expression by truncating the operators Fi tk (N) to FL tk (N ) as 
described above, then: 

Z kj^O zv z 

Write out: 

W v ( k ) F LANo)F L ,.k(N ) 

Z k^O 

= 7W7 X v ( k ) ( a L,o a l,o a L,k a L,k + a* Lfi a Lt0 a Lt - k a* L _ k 
zv k^O 

+ a L,0^L,0 a *L t k a *L-k + a *Lfl a *L,0 a L-kO'L,k 

= W77 X v ( k ) (( a L,oa*L,o + a l,o a L,o)a* Ltk a Lik + a Lfi a Li0 a* L 

ZV k^O 

j_ * * \ j_ /v 

+ a L,O a Lfi a L-kaL,k) H — iv L)0 

a" 

As in [23|, replace the operators by complex numbers |c|e =l=, ' a: in this part 



k a L,-k 



of the interaction and preserve them as operators in the remaining terms, 
then 

Ul = \Y,<m\e- 2ia <,k<,-k + z 2ia ^-ka L ,k) 

+ M 2 X <k)at tk a Lik + V M N 2 _ l mNh + ffi, c | V. (81) 
fc^O zv z z 

Apart from the term |0(O)A^, which only leads to a shift in the chemical 

Ho 

2 



potential, and an unimportant constant ^^-|c| 2 V this is exactly the Hamil- 



tonian (p5f). 



4.2 Collective Goldstone modes 

We consider again our system to be in one of the extremal equilibrium states 
ujg, and without loss of generality we take a = (i.e. c real), and denote 
this state by ujp. As before, let 

F L , q {N ) = (a* L q a Lfl + a* L a L ^ q 

F L , q (A) = -^ 2 jy( x )-a{x))e^dx 

= -^V a L,q - a L,-q), 
21 



and for ease of notation: 

P° L , q = 7p F U^o), A L , q = F L , q (A). 
These operators still satisfy the correct commutation relation 
\pl, q , A l ,_,] = ^j== (a* L>0 + a L ,„), 

and by ([7|): 

hm{pl g ,A L ^]=i. (82) 

More generally 

Um[^,,A w ]=iV- (83) 

i>— >oo ^ 

Proposition 8 We /iai>e /or q,q' 
(i) 



^,/w) = V^-coth(^ 



hm u p {p Uq p L _ ql ) = 8 q , q > — coth 

lim = ^-coth — 2 . 

L ^°° 2e g \ 2 J 

Proof. This is an easy calculation using the Bogoliubov transformation ( 
(|68|), the explicit form of the state ([Q]), property (|75|) and the fact that 



cosh 2a q 



sinh 2a q 



e q + c 2 v(q) 



E q 



c 2 v(q) 



so that 



(cosh a q + sinh a q ) 2 = cosh 2a q + sinh 2a, 



e, + 2cM<?) 



54) 



22 



(cosh a q — sinh a q ) 2 = cosh 2a q — sinh 2a q 

e q + 2c 2 v(q) 
E q 

El 

e <? 

le g + 2c 2 v(q) 

= v~ — • 

This then gives the result via 

= &q,q' t: (cosh a q + sinh a q ) 2 X 



hm up ((b* Lq + b L ^ q )(b* L ^ q + b L , q )), 




□ 

For p < oo the small g-behaviour of these variances is 

lim uj p (p° L p\ ) ~ const ^ 
lim Uf3(A Lt gA L - g ) ~ const x |g|~ 2 . 

So again we have the phenomenon that at non-zero temperature it is impos- 
sible to do a renormalisation of p° L and ? which gives both a meaningful 
g — > limit for the variances and preserves a non-trivial commutation rela- 
tion. 

Therefore we will restrict ourself from now on to the ground state (denoted 
a;). Contrary to the imperfect Bose gas, there remains a non-trivial in- 
dependence in the ground state. This is because even at zero temperature 
not all particles condense into the ground state. We have for small q 

lim u(p° L p° L ) oc \q\ 
lim up(A Lq A L q ) oc 

L— >oo 

Remark that for the condensate density fluctuations this is the above men- 
tioned linear behaviour. 
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We now redefine p\ and A^q to be self-adjoint and renormalised in q, i.e. 



PL 



l L,q + a L,~q) a L,0 + a*L, ( a L,q + ^L,- 



A 



2{c 2 \q\V) 1 ' 2 L v 

|l/2 

a L, q + a L,- q - ( a L,q + a L - q ) 



\Q\ 



(86) 
(87) 



Hence we have still 



lim \PL,q, A L,q>] = iS, 



L— >oo 



q,q ■ 



It is already interesting at this stage to remark that 



lim lim u ({PlJ 2 



lim lim u p [A L 

q— >0 L-+00 v 



1 

2ft 



with 



And hence 



Vt = lim 



■ E q \q\ 



(4mc 2 u(0)) 1/2 . 



n 2 lim lim u ((plj 2 ) = lim lim u (A 2 L 



(89) 
(90) 

(91) 
(92) 



Again the rest of this section is devoted to the rigorous mathematical treat- 
ment of the existence of the fluctuation operators. Again at first reading, the 
reader can immediately proceed to section fO| 

With all definitions and notations as above for the imperfect Bose gas, we 
define again for f,g G T 



Plq(f) 



/(?> o)01 q + a*L, -qWfi 



and 



2cV x / 2 

+/(0,g)a2 (a Lig + a L> 

% r 

9(q, 0) (aL,q + a *L,-q) ~ 9(0, q) {a L>q + a L> _ 9 ) 



F L>q (f,9)=pL(f)+ A LM- 



(93) 
(94) 



(95) 



Now we take immediately the double limit lim 9 _>o limL->oo rather than the 
two limits separately as for the imperfect Bose gas. 
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Proposition 9 For f^gEj 7 , 

lim lim u {F Lt q(f, : 



lim 



e q + c 2 v(q) 



+o 2E„ 



f(q,0)+ig(q,0) 



2 c 2 v(q) 
~ 2E n 



5R 



(/M) + ^(g,o) 



(96) 



Proof. An explicit calculation shows that 



lim u> ( P i q (fy 



lim a; (A L>q (g)' 



e g + c 2 t> (g) 
2£ g 

e g + c 2 v{q) 
2E n 



g(q,o)\ 2 + 



c 2 v(q) 

c 2 w(q r ) 
4£„ 



[/M) 2 + /(0,g) ; 



[g{q, 0) 2 + g(0,qf 



and 



f ' 1 r{q) -W(q, 0)g(0, q)] + ^9[f(g, 0)], 



E n 



which together give (|96[) . 



□ 



From now on we restrict to those pairs of functions (/, g) for which (|96|) 
remains finite. This means if / is real, \f(q, 0)| should not diverge faster 
than |g| _1//2 and if / is imaginary, 1/(5,0)1 should go to zero, at least like 
Igl 1 / 2 . A real g should satisfy the same condition as an imaginary / and vice 
versa. Alternatively, we could say we restrict T to those functions / which 
satisfy the above conditions and then look at pairs (/, Jg), with /, g in (the 
restricted) T and Jg defined as (Jg)(q,0) = —ig(q,0), (Jg)(0,q) = ig(0,q). 



Proposition 10 (A BCH-formula) Let (fi,gi) G (J r ,JJ-'),i = 1,2, then 



lim 

L— >oo 



JF L , q (h,gi) e iF L , q (fr>92) _ e i{F L , q {h^i)+FL, q (hm)) e -\[FL, q {hm),F L , q U^92)] 



0. 



(97) 
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Proof. See Appendix [A]. 



□ 



On V, the space of complex linear combinations of elements from (JF, J J 7 ), 
JF in the above defined sense, define a positive sesquilinear form < -|- > by 

< hi9i\h,92 >= lim lim u(F Ljq (f l ,g i yF L , q (f 2 ,g 2 )), 

and extension to V by linearity. This can be calculated in the same way as 
(|96|) is calculated to give the following expression 

lim u(F Lt q(f 1 ,g 1 )F Liq {f 2 ,g-i)) 



e q + c 2 v(q) ^ fyj : — Trj. . ,\ c 2 v(q) 



2E q • " "' ^ '2 /-.',, 



+ i$([f l + i gi }[f 2 + ig 2 }). 



where we adopted the short hand notation /j = fi(q, 0), fi = fi(0, q), etc. 

In practice all operators we use, have either / = or g = 0, and / and g 
will always be either real or imaginary. It is easily seen that this leads to a 
considerable simplification of the above formula's. 

Separating the real and the imaginary part of the restriction of < -|- > to 
the real subspace (F, JF) of V, 

< fx,9\\hi92 >= 9i\f2,92) + ^(fi,gi\f2,92), (98) 

defines a real bilinear positive symmetric form s and a symplectic form a. 

Remark that, as it should, this a is the same as the one for the imperfect 
Bose gas since it also satisfies 

lim lim [F Liq (f l7 g t ), F L:q (f 2 , g 2 )\ = ia(ft, gi\f 2 , g 2 ), (99) 

q— »U L— >oo 

where the limit is taken in the GNS-representation of the equilibrium state 
uj. This commutator is of course independent of the given model since in 
both models we have 

V- l/2 ai - c. 
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Proposition 11 (Central Limit Theorem) For (f,g) G (J 7 , J J 7 ), t G R 
lim lim w fe^.^/.f)) = e -£'(f*\M. ( 10 o) 

Proof. Using (|75|) we get 

lim lim u ( e itFL ^ f ^ 



lim lim u (e% {lf(*>W+ i 9(<ifi)]( a t , q +<,- q )+lf(^)+ig(qM(aL, 9 +a L ,- q )}\ 

With the explicit expression of the state (|7l|), the r.h.s. of (p. 01] ) can easily 
be computed and we get ( |100| ). 

□ 



Proposition 12 (Reconstruction Theorem) The linear functional 

Q(W(f,g)) = e ^(/> 9 |/, 9 ) ( 102 ) 

defined on the algebra CCR^J 7 , J J 7 ), a), is a quasi free state. 
More explicitly, we have for all (/,, gi) G (J 7 , J J 7 ), i — 1, . . . , n, 



lim 



9- 



i lim w f e ^.«(/i.ff0 . . . e ^Afn,9n)\ = q iWlf^gx) ■ ■ ■ W(f n , g n )) . (103) 

L— >oo V / 



The state u is regular and hence for every (f, g) there exists a self-adjoint 
Bo sonic field $(f,g) in the GNS representation 7T£, f^) such that 

*o(W(f,g)) = e i *V>«\ ( io4) 

This implies that in this sense of the central limit ( \1 0$) , the local fluctuations 
converge to the Bosonic fields associated with CCR^J 7 , J J 7 ), a): 

$(f,g) = CLT - lim lim F L , q (f,g). (105) 

q^O L~^oo 

Proof. See the proof of Proposition | and Proposition |6|. 

□ 
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Specifying to our original operators: 



P° = $(|g|- 1/2 ,0) = CLT-limlimp^ 

q^O L— +00 

i = $(0, |g| 1/2 ) = CLT-lim lim A Lq , 

q^O L— >oo ' 

adopting the notation (f,g) = (f(q, 0),g(q,0)). 

The equivalence property also remains true in this case of course, i.e. 



(106) 
(107) 



®(fi,9i) = $(/2,#2) 



(/l,#l) ~ (f2,92) 

lim lim w (F £) ,(/i - / 2 , #1 - # 2 ) 2 ) = 0, 



in particular 



CLT - lim lim pl q (f) = CLT - lim lim A Lq (Jf). 



q^>0 L— >oo 



5^0 L— »oo 



(108) 



4.3 Dynamics of the collective Goldstone modes 

We will now show that also for the weakly interacting Bose gas the fluc- 
tuations of the generator and of the order parameter of the SSB decouple 
dynamically from the other degrees of freedom of the system, and form a 
harmonic oscillator system. 

At the local level, we have 



2(c 2 \q\Vy/ 2 
^ ic 2 v(q) 



Pi, 



2(c 2 \q\Vy/ 2 



[<, q + a* L ,-q) a L,o ~ a* Lfi (a Liq + a L ^ q ) 

[ a L,q + a *L,- q )( a L,0 - 4,o) + K,0 - a*L, )( a L,q + dL-q) 

ic 2 v(q) 



+ 



ql 1 / 2 ) 2(c 2 \q\Vy/ 2 
+ (a L ,o - a* Lfi )(a L , q + a L ^ q ) 



(109) 



Since the operators V 1//2 (ol,o ~~ a *Lo) converge to in the GNS represen- 
tation of oj, we only have to take into account the term p\ q ( 7^172 )■ This 
is what we mentioned earlier, the interaction part of Hl(c) commutes with 
the condensate density fluctuations as one expects physically, but only in the 
thermodynamic limit. 
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Also: 



-i[H L (c),A L>q ) 
(e g + 2c 2 v(q))\q\ 1/2 



a L,g + a L,-q + a L,g + a L, 



A L , g (-i(e g + 2c 2 v(q))\q\ 1 



/2 



(110) 



As in the imperfect Bose gas, we define a macroscopic dynamics by first 
renormalizing the Hamiltonian with its spectrum, and then taking the limit. 
CLT — lim^o hniL-^oo of the commutators. 



i[H,p 



~0l 



CLT — lim lim i 

q— >0 L— >oo 



TrH L (c),pl 



CLT — lim lim p ° r „ , 

9-0L^oo FZ "9 \E q \q\ 1 / 2 1 

CLT — lim lim A^ A \ — j t \q 



1/2 



q— >0 L— >oo 



(in) 



(where we used the equivalence relation). Remember that Al, 9 = Ai^Qql 1 ' 2 ) 
and that Q was defined in (|9T|). 

n = (4mc 2 t;(0)) 1/2 . 



- zfF, Al = CLT - lim lim -i 

1 1 q^OL^oo 



-H L {c),A Liq 



CLT — lim lim Al„ 

q— >0L^oo 



-i{e q + 2c 2 v{q))\q\ 1 l T 



E n 



Eq 


q\ 


i 




\q 


1/2 



= np u . 

Again remember p\ = p^ Mq]" 1 ^ 2 )- 

Hence we have found a canonical pair p° and A, satisfying 



112} 



;ii3) 
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which dynamically decouple from the other degrees of freedom of the system, 
with the dynamics given by 

i[H,?\ = ~A (114) 
-i[H,A} = np°. (115) 
The solution of these equations is the harmonic oscillator with energy Q: 

H=±(tf{jP? + A*). (116) 

The state u is an equilibrium ground state for the Hamiltonian H satisfying 
the following virial theorem 

Q 2 u) ((p ) 2 ) = u (I 2 ) . (117) 



We classified the programme of Anderson JT3|, [TJ| about the construction of 
the canonical Goldstone coordinates. 

Remark that also in this case we have done a rescaling of time 

t 



t -> t 



E q 



when going from the microdynamics Hl(c) to the macrodynamics H, indi- 
cating again that the typical lifetime of the fluctuations becomes infinite as 
q — ► 0, this time with a rate \q\~~ 1 , as a consequence of the fact that E q oc \q\ 
for small q. 
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Appendices 

A Proof of Proposition [3] and |tO| 

In this proof we make use of the following general formula 

e it(x+y) = e itx + dse isx iye i{t - s){x+y) , (118) 
Jo 
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where x, y are self-adjoint operators. From this formula, one deduces useful 
properties like for example 

[e ix , z]=i f 1 dte itx [x, z}e i{1 - t)x , (119) 



\\K,z]\\<\\[x,z]\\, (120) 
where x,y are self-adjoint, z arbitrary and |.| is some norm. 
Denote = Fjj ig (fi, gi), i = 1, 2. We need to prove 



One has 



lim \\e m e iF * - e *(*i+*» e -*[*i.*y || w = 0. 

L^oo 



\\e e — e v 'e 



^| e iFi e iF 2 _ e i(F 1 +F 2 ) e -^[F u F 2 ]y^ e iF le iF 2 _ e i(F 1 +F 2 ) e -^[F 1 ,F 2 \y 

= 2-lu (e- iF2 e- iFl e< Fl+F ^e-^ Fl ' F ^) - cu (e^ F ^ e - i{ - FM) e iF 'e iF2 ) . 
Hence it is sufficient to show 

lim u (e- lF2 e- m e l{Fl+F2) e-^ F2] ) = 1. 

Define a function f(t) by 

f{t) = u (e-^e- im e^ +F ^e-^ F ^) - 1. 
We have to show 

lim |/(1)| = 0. 

By Taylor's theorem, there exists some t e [0, 1] such that 

f(l) = f(0) + f(t) = f(t), 
since /(0) = 0. Let us calculate f'(t) using ( |118| ), 

J(tF 1 +F 2 ) = i f 1 ds( ,is(tF 1 +F 2 ) F ^ e i(l- s )(tFx+F 2 ) _ 
JO 



d 

— e" 



Hence 

+ie~ iF2 e~ itFl rf se is( * Fl+F2) Fie i(1 ~ s)( ' Fl+F2) e~^ [Fl ' F2] 
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Denote 



[Fi,F 2 ]e 

+2e- iF2 e- i * Fl F 1 e i(tFl+F2) e-^ [Fl ' F2] 

+ie -iF 2e -it Fl f 1 ds { e is(tF 1+ F 2 )^ F ^ e i { l- s ){tF 1+ F,) e -l[F u F 2 ] 

Jo 

_ I e -^2 g-itFx +F 2 ) ^ 2 ] e -| [*i ,F 2 ] j 

y 1 ds {[ e -(^i+^) )Fi ] e i(l-,)(tF 1+ F 2 ) 

A = £ ds |[ e -(^+^), Fl ] e *(i-s)(tF 1+ F 2 ) - ^ Fl+F2 ) [F 2 , Fx]} . 



= u [ie-^e-^e-^'^A + ie'^e'^ [A, e~f ^ 



then 



|/'(*)| < |w (e- lF2 e-^e-5^^]A)| + \u (f iF * e' itF ' [A, e'^^ 
The first part is estimated as follows: 

^(e-^e^e-s^U)! < ||A|| W , 
by the Cauchy-Schwarz inequality. We now make an estimation of 

£ ds {[ e is ^ +F *\ Fl ] e ^-^+F 2 ) _ ise i{tF,+F 2)[F ^ Fi] j 

< /"'ds [ e »^+^) ) F 1 ] e 'M(^+^)_^+^)[ F2)F i 
Jo llw 

= /'ds /'dre^^+^^stFa.Fxle^-^^+^Je^ 1 -')^ 1 ^) 
jo jo 



(121) 



PL = 



-ise 



= j 1 ds f 1 dre irsi - tFM \s[F 2 , Fi y(i-rs)(tF l+ F 2 ) 

JO JO 

-ise i ^ +F ^[F 2 ,F 1 ] 
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< ds drs 
Jo Jo 



r/.s / drise irs ^ +F ^[[F 2 ,F 1 },e i ^- rs ^ +F ^ 

[F 2 ,F 1 ],e i{1 - rs){tF ' +F2) 

< f 1 ds f 1 drs{l-rs)\\[[F 2 ,F 1 ],tF 1 + F 2 
Jo Jo 



^\\[[F2,F 1 ],tF 1 + F 2 ]\\ ui 



where we have used (|119|) in the third step and ( |120| ) in the fifth. But because 
the commutator [F 2 , Fi] converges to a complex number in the state u, it can 
be calculated explicitly that the commutator [[F 2 , F 1 ],tF 1 + F 2 ] converges to 
0. By a similar argument it can also be shown that the second part of (|121|) 
converges to 0, thus proving Proposition [3[ 



□ 



B Proof of Proposition [4| 

Take some arbitrary FL >q (f,g) and first of all remark that it can be written 



as 



F L>q (f,g) = B; + B q , 

where B* = B^ q and [B*,B q ] = (we have only written the g-dependence 
explicitly, all other dependencies are implicit). One can for instance take 



B„ 



1 



2(p V) 1 P 



22 f(k + q, k)a* L>k+q a Lik + - [g(q, 0)a* L q - g(0, q)a L> 



B„ itself can be decomposed into 



B q = B q + B q , 



:i22) 



with 



1 / 



f{q,0) 



O'Lfl 



(Povy/i 



ig(q,0) a 



L,q 



+ v (0 ' g) (p vy/ 2 ~ g) j aL ' 
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We want to show now that the part {B q )* + B q is the only part of g) 
which gives a non-zero contribution to the expectation value on the l.h.s. 
of (P^). Expanding the exponential in a power series, we have to calculate 



expectation values of the type 

{(b; + B q 



This is obviously zero for m odd, and for m = 2n even, the only non-zero 
terms are those with a number of starred operators equal to unstarred. Be- 
cause of the commutation [B*, B q ], these are all equal to 

u> ((B*) n Bfj . 

Using the decomposition ( |122|) this becomes 

u ({B* q ) n B r q l ) = u ({{B° q )*) n {B° q ) n ) + other terms. 

We now prove that these 'other terms' are all necessarily zero because we are 
working in the ground state. Using the commutation of the (-8°)" with the 
{B q )* these terms are all of the form 

" ((^ )*) i (^ ) i (^) n " i (^)"" i ) • (123) 

We have 

(B q T~ j = - E f( k i + 9, *i) ■ ■ ■ f(kn-j + Q, kvr-i) 

(4p U) 2 klr .. >kn _..- q ^ kl ^o 

Using this in ( |123|) one gets sums of expectation values which can be com- 
puted by using the quasi-freeness of the state u, i.e. each expectation is a 
sum of products of one- and two-point correlations, respectivily ui{a Lk ) and 

of which only the following are non-zero: 



^(4,o) = \[p<y 

uF{a Ltk a* L k ) = u{a Ltk a* L k ) = 1, k ^ 0. 

We show now that each of the expectations to be computed is zero. Take an 
arbitrary term. First of all, take some ai,<aa* Lq term from one of the B q 's. It 
can only give a non-zero contribution if it is combined with an a* L cil,<j term 
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from one of the (-B°)*'s, so at least we need i = j. A term ± from any of the 
.Bg's or (i?°)*'s can never be combined with a term coming from the (B q )^s 
because this would give rise to an expectation of an odd number of creation 
and annihilation operators, all with a non-zero index, which is zero. Hence in 
the quasi-free decomposition, the operators coming from ((B®)*y (B®y and 
the operators coming from (B*) n ~ : > (5 (? ) T1_ ' J completely decouple from each 
other. But a typical factor arising from (B*) n ~ : > (B q ) n ~ : ' is 



uo a 



l L,k 1 +q CLL M - - ' a L,k n -j+q a L,k n _ j a L,l 1 -q a L,l 1 ■ ■ ■ a L,l n _j+q a L,l n „ 

which is zero because it always contains at least one factor uj(a* L ^ai^), k ^ 0. 
Hence we proved that indeed 



uj ( (B* q ) n B n q 



or in other words 



hm u U tF ^9)\ = i im u ( e *n«U>ffi) , 



(124) 



where 
with 

P°L, q (f) 

Hence 



FlJf,g) = p° L Jf) + A L M, 



/(?>0) / * . * \ , /(o,g) * f 

Yj^ a L,o [a L q + a L _ q J + J—yyJ5 a L,0 K a L,q + dL,- 



(povy/ 



F L , q (f,g) 



2 ' 

+ 



f(q,0) T ^ T7 , + ig(q,0) 



(Povy/ 2 



f(P,q) 



b L,0 



(p vy/ 2 



{0-L,g + a L,-q 
0>L,q + a L-q) 



Using ( |T9| ) one gets 



lim uj [ 



e itF °L, q (f>9) 



lim uj ^{ [/(9 ' 0)+ * 9{9,0) K a ^ +a £,-J^ . (125) 

With the explicit expression of the state fll6|) , the r.h.s. of ( |125| ) can easily 
be computed, and together with (|124|) one gets (|39f) . 



□ 
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C Proof of Proposition [5| 



We prove (f4T|) by induction onti G No- The basis of the induction is the 
Central Limit Theorem, Proposition f|. 



To prove the induction step, assume that holds for some n G No and fix 
(fi,9i) e (F,F),i = l,...,n + l. 

For convenience write 



e iFL,q(fl,gi) . . . e iFL,q(fn-l,gn-l) = ffi 



and 



FL, q (fh9i) = F l Lq . 
By the Cauchy-Schwarz inequality and the BCH-formula 



lim 

L^oo 

< lim 

L^oo 



10 



(w L , Q 



' iFV- iF" +1 i(F? +FV-+ 1 ) 

p L,qp L,q — p V L,q 1 L,q >p 2 L ' 



iF r n iF r n+1 i(F? +FV-+ 1 ) -h\F? ,F? +1 )] 

p h,qp L,q — p ^ L,q 1 L,q 'p 2 L L,q 'J 



Use the Cauchy-Schwarz inequality again to derive that 



"'L.+^e-sI^,,.^ 1 )]^ _ 



L(?e '( F r,«+ F r,« 1 )) e -i^Hn+i) 

2 



_ 1 [PTl + 1 

g 2 y r L,q-> r L,q 



2 -w ( e ^ [F ^' F ^ Jl )e"5 



jCTq(n|n+l) 



e -50-,(n|n+l)j 



This expression converges to zero as L — ► oo because of the convergence of 
the commutator. Combining the induction hypothesis and the above result, 
one finds that 



\imu(W Lq e iF ^e iF ^) = Hmu(w Lq e^- 

(W 9 (l).--[W^(7i) + W 9 (n + l)] 



e 2 



§<r a (ra|n+l) 



= w«(W,(l)...W,(7i)W,(n + l)). 

The last equality results from the CCR algebraic structure of 
CCR((F,f),<r q ). 
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The only thing left to prove is positivity. Take f, gi 6 J 7 , i — 1, 2 and use the 
definitions of a q and s q along with the Cauchy-Schwarz inequality to derive 
that 

4 W q (fl,9l\f2,92)\ 2 < \< fl,9l\f2,92 >q\ 2 
<< fl,9l\fl,9l >q< f2,92\f2,92 >q < $q (fl , 9l I fl , 9l) $q (/a, 92 \ jfe, 92) ■ 

□ 

D Proof of Proposition |7| 

Denote 

W(f i ,g i )=W i . 
Suppose first that (ii) is satisfied, then 

[^(W 1 )^(W 2 )\ =0 

and hence 

cr(l|2) = 0. 

Further 

1 = $ (W X W%) = Q (W X W„ 2 ) 

where we used the notation (fx — f 2 , g\ — g 2 ) — > 1 — 2. By the definition of s 
this means < 1 — 2|1 — 2 >= 0, proving (i). 

Conversely, suppose (f x ,gi) ~ (f 2 ,g 2 ) then 

~ \a(l - 2\x)\ 2 << 1 - 2 1 1 - 2 >< x|x > 

implies that cr(l — 2\x) = for all x, where x denotes an arbitrary element of 
(J 7 , J 7 ), i.e. tt^ (Wx- 2 ) commutes with all elements of n^(CC 'R((J- ', -T 7 ), er)) = 
A4, or ttq (W1-2) belongs to the commutant A4' of M.. Also 

\\(*g,(w x - 2 ) -i)^f = w^-in^-i)) 

= 2-w((Wi_ 2 )-a;((W 2 _i) 
= 0. 



37 



As is cyclic for M. it is separating for M' . Hence 

TTu, (^i_ 2 ) = 1 

or 

n a (W 1 )=n S) (W 2 ). 

□ 
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